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J^j- abstract 

r j Aims. Investigation of the dense gas, the outflows and the continuum emission from the massive twin cores NGC6334I and I(N) at high spatial 

^ ' resolution. 

0O Methods. We imaged the region with the Australia Telescope Compact Array (ATCA) at 3.4 mm wavelength in continuum as well as 
C — CH 3 CN(5k - 4 K ) and HCN(l-O) spectral line emission. 

Results. While the continuum emission in NGC6334I mainly traces the UCHii region, toward NGC6334I(N) we detect line emission from 
, four of the previously identified dust continuum condensations that are of protostellar or pre-stellar nature. The CHaCN(5x - 4 K ) lines are 
y—{ • detected in all A"-components up to energies of 128 K above ground toward two protostellar condensations in both regions. We find line- width 
(*^") increasing with increasing A* for all sources, which indicates a higher degree of internal motions of the hotter gas probed by these high 
OO , K-transitions. Toward the main mm and CH 3 CN source in NGC6334I we identify a velocity gradient approximately perpendicular to the 
■ large-scale molecular outflow. This may be interpreted as a signature of an accretion disk, although other scenarios, e.g., an unresolved double 
source, could produce a similar signature as well. No comparable signature is found toward any of the other sources. HCN does not trace the 
• i— I . dense gas well in this region but it is dominated by the molecular outflows. While the outflow in NGC6334I exhibits a normal Hubble-law like 
' velocity structure, the data are consistent with a precessing outflow close to the plane of the sky for NGC6334I(N). Furthermore, we observe a 
t wide (~15.4km s -1 ) HCN absorption line, much broader than the previously observed CH 3 OH and NH 3 absorption lines. Several explanations 
. . . i for the difference are discussed. 
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1. Introduction 

The massive twin cores NGC6334I and I(N) 



at a 



dis- 



tance of 1 .7 kpc in th e southern hemisphere (Neckel 1978; 
Straw & Hvlandlll989l) have been subjected to investigations 
for more than two decades. The two regions are located at the 
north-eastern end of the much larger molecular cloud/Hn re- 
gion complex NGC6334 ( e.g. 



1982t IGezari 

1982; de Pree et all 1 19951; iKraemer & Jackson! 1 19991; ISandelll 



Rodriguez et al 



200rilCarral et alJl2002l) .The reason why NGC6334I (synony 
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mous with NGC6334F) and NGC6334I(N) are so interesting 
from a comparison point of view is that they are only sep- 
arated by approximately 1 parsec, hence they share a similar 
large-scale molecular environment, but they exhibit extremely 
different characteristics likely because they are at different evo- 
lutionary stages. 

Both regions have been studied in much detail over the 
last decades; rece nt summarie s of th e past observat i ons can 



be found, 



Ml: 



Rodriguez et al. 



Junter etalJ d2006h . iBeuther et alJ J2007bl) 
d2007t) . Here we just outline their main 



characteristics. NGC6334I is a prototypical hot molecular core 
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right at the head of a cometary ultracompact Hn ( UCHii) 
region ddePree et alJll995t: iKraemer & Jacksonlll995|) : It ex- 



Table 1. Line parameters 



hibits rich spectral line emission jMcCutcheon et al. 2000; 



IThorwirth et al.l 120031: ISchilke et al.ll2006l). a bipo lar outflow 
dBachiller & Cernicharolll990l: iLeurini et alj|2006h and H 2 0, 
OH, CH3OH class n and NH^(3,3)/(6,6)/(8,6)/(ll,9) maser 



emission (Moran 



Gaume & Mutel 

"Tr 



& Rodriguez 1980; Forster & Caswell 1989; 
Brooks & Whiteoaldl200lUNorris et al l 



198_p 



1993; Caswell ! 1 19971: iWalsh et alJU998t iBeuther et alJl2007bt 
Walsh et alj |2007). In contrast to that, up to very re- 
cently NGC6334I(N) was considered a typical cold core 
since no mid-infr ared and only faint near-infrared emis- 
sion w as detected (lGezarilll982t iTapia et a l. 1996; Persi et alJ 
2005). Furthermore, weak cm continuu m and class i and 
ii CH3OH maser em i ssion was reported dCarral et alj 120021: 
Kogan & Slvshlll998l; ICaswelll 1 1 997t IWalsh et al.lli998l) . The 



spectral line forest is considerabl y less dense compared to 
NGC6334I (IThorwirth et alJ 120031) , however, a few species 



are stronger toward NGC6334K N) dSollins & MegeathlH oa 
Walsh et al. in prep.). In addition. iMegeafh & Tieftrunkl(ll99' 
report the detection of a molecular outflow in this region as 
well. In summary, both regions show signs of active star for- 
mation, however, the southern region NGC6334I appears to be 
in a more advanced evolutionary stage than the northern region 
NGC6334I(N). 

To better characterize this intriguing pair of massive star- 
forming regions, we started a concerted campaign from cm 
to mm wavelengths with the Australia Telescope Compact 
Array (ATCA), the Submillimeter Array (SMA) and the Mopra 
single-dish telescope. The previous ATCA NH3(1,1) to (6,6) 
line observati ons revealed compact war m gas emission from 
both regions ( Beuther et al.l2005 , 2007bl) . and temperatures es- 
timated to exceed 100K. While toward NGC6334I(N) the low 
energy NH3 lines showed only extended emission, the high 
energy lines finally revealed compact gas components. The 
NH 3 (6,6) line profile from NGC6334I(N) allowed speculation 
about a potential accretion disk. CH3OH was strong in absorp- 
tion toward the southern UCHn region in NGC633 4I, indicative 



of exp anding gas. In the mm continuum emission, Hun ter et al 



(2006) used the SMA to resolve several mm continuum sources 
toward both regions (4 in NGC6334I and 7 in NGC6334I(N)). 
Furthermore, Hunter et al. (in prep.) identified an additional 
SiO outflow in NGC6334I(N) that has its orientation in north- 
east south-west direction, approximately perpendicular to the 
one previously reported bv lMegeath & Tieftrunk ( 1999). 

Here we present 3.4mm continuum and HCN/CH3CN 
spectral line observations obtained with the new 3 mm facil- 
ity at the ATCA. These observations shed light on the outflow 
and dense gas properties of both regions as well as on the con- 
tinuum emission from the embedded protostellar objects and 
the UCHn region. 

2. Observations 

The two regions NGC6334I & I(N) were observed in May 
2006 during two nights with the ATCA in the H214 con- 
figuration that results in projected baselines between 14 and 
73 kA at 88 GHz. The weather conditions at Narrabri were ex- 
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[10 5 cirr 3 ] 


HCN(l-O), F=l-1 
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HCN(l-O), F=2-l 


88.632 


4 


32 


HCN(l-O), F=0-1 


88.634 
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32 


CH 3 CN(5 4 - 4 4 ) 


91.959 


128 


7 


CH 3 CN(5 3 - 4 3 ) 


91.971 


78 


7 


CH 3 CN(5 2 - 4 2 ) 


91.980 


42 


7 


CH 3 CN(5! -40 


91.985 


20 


6 


CH 3 CN(5 - 4o) 


91.987 


13 
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For all observed lines the columns list (left to right) the 
species/quantum numbers, frequencies, upper level energy states E u /k 
(k is the Boltzmann constant) and critical densities n crit . 
" The critical densities n cl j t = A/y are calculated at 60 K (Einstein co- 
efficient A and collisional rate y ). For HCN, A and 7 are taken from 
LAMBDA dSchoieret all 120051) . For CH 3 CN, A is calculated from 
A = 0.3/l^ yu 2 (/lioo in units of 100 /mi and n = 3.9 debye), and the 
corresponding y-values are from lPei & Zend jl995bl [a| ). For a line op - 
tical depth r > 1, ra cr j t has to be multiplied by 1/t (Tielens 2005). 



cellent with approximate precipitable water vapor of -10 mm 
and measured system temperatures between 150 and 670 K. 
The phase reference centers were R.A. (J2000) 17' 1 20" , 53 s .44, 
Deck (J2000) -35°47'02".2 for NGC6334I and R.A. (J2000) 
17''20 m 54 s .63, Deck (J2000) -35°45'08".9 for NGC6334I(N). 
The velocities relative to the local standard of rest (vj sr ) for 
NGC6334I and NGC6334I(N) are ~ -7.6 and ~ -3.3 km s 1 , 
respectively. We observed the 3.4mm continuum emission at 
88.4 GHz with a bandwidth and spectral resolution of 128 and 
1 MHz, respectively. The spectral range for the continuum was 
checked to be line-free based on previous MOPRA single-dish 
observations of that region (Walsh et al., in prep.). During one 
night, simultaneously with the continuum emission, we ob- 
served the HCN(l-O) line at 88.632GHz. In the second night, 
the CH 3 CN(5a: - 4 K ) transitions at -91.98 GHz were targeted 
averaging two polarizations to achieve better signal-to-noise ra- 
tio. For more details on the spectral lines, see Table [T] A good 
uv-coverage was obtained through regular switching between 
both sources and the gain calibrators 1742-289 and 1759-39. 
The channel separation of the HCN observations was 0.25 MHz 
(-0.85 km s _1 ) and slightly worse with 0.5 MHz (-1.63 km s -1 ) 
for CH3CN because we had to cover a broader bandwidth due 
to the several -components. Flux calibration was performed 
with observations of Uranus and is estimated to be accurate 
within 20%. The primary beam of the ATCA at the observing 
frequency is 36" (FWHM). The data were reduced with the 
MIRIAD package. Applying different weightings for the line 
and continuum data - mostly uniform weighting for the com- 
pact continuum and CH3CN emission and natural for the more 
extended HCN emission - the synthesized beams vary between 
the different maps. The achieved spatial resolution and lcr rms 
values are given in Table [2] 
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Table 2. Synthesized beams 9 and rms 



Line 


G 


lcr rms 




["] 


r mJy -i 
beam ^ 


NGC6334I 


3.4mm cont. 


2.3" x 1.6" 


23 


CH 3 CN(5 4 - 4 4 ), [-12,2]km/s a 


2.2" x 1.6" 


7 


CH 3 CN(5 4 - 4 4 ), 1.7km/s i 


2.2" x 1.6" 


11 


HCN(l-O), [-25,-19]km/s a 


2.6" x 1.8" 


19 


HCN(l-O), [5,19]km/s a 


2.6" x 1.8" 


16 


HCN(l-O), 1.7km/s* 


2.5" x 1.8" 


33 


NGC6334I(N) 


3.4mm cont. 


2.6" x 1.8" 


2.6 


CH 3 CN(5, =al -4 t=0A y, [-12,8]km/s 


2.2" x 1.6" 


4.9 


CH 3 CN(5 4 -4 4 ), 1.7km/s* 


2.2" x 1.6" 


8.4 


HCN(l-O), [-16,-6]km/s a 


2.5" x 1.8" 


19 


HCN(l-O), [4,12]km/s fl 


2.5" x 1.8" 


15 


HCN(l-O), 2.0km/s i ' 


2.5" x 1.8" 


26 



" Integrated velocity regime. 

* Channel maps with the given velocity resolution. 



3. Results and Discussion 

All of the targetted spectral lines with excitation levels above 
ground, E u /k, between 4 and 128 K (Table [TJ and the 3.4mm 
continuum emission were detected and mapped toward both 
target regions (Figs. |T|&[2ji. While CH 3 CN is a typical hot core 
molecule and only detected toward the central warm protostars, 
HCN exhibits significantly more extended emission and shows 
a strong association with the various molecular outflows in the 
two regions. This may be considered surprising since the crit- 
ical density of HCN(l-O) is even larger than that of CH3CN. 
The 3.4 mm continuum emission is of different origin in both 
regions: while it is largely due to the free-free emission from 
the UCH11 region in NGC6334I, toward the northern region 
NGC6334I(N) the continuum emission stems mainly from the 
dust in the vicinity of the embedded protostars. In the following 
we will outline the characteristics of both regions separately. 

3.1. NGC6334I 

3.1.1. 3.4 mm continuum emission 

Figure Q] presents an overlay outlining the main features of 
the spectral line and continuum emission toward NGC6334I. 
Almost all of the 3.4 mm continuum emission arises from the 
UCHn region with a comparable morphol ogy to the previous 



cm continuum images of that region (e.g.. lde Pree et al. 1995; 



Beut her et al.ll2005l) . In contrast to that, it is only barely de- 
tectable at a 3.8cr level of 87mJybeairr' toward the strongest 
1.4 mm continuum, NH 3 and CH3CN peak positio n, mml 



in Fig. Q] dHunter et alJ l200fjt iBeuther et aD |2007b), mark 



ing the location of the dominating protostar in the region. 
Comparing this 3.8< x detection with the 1 .4 mm data-point from 



Hun ter et alJ (J2006) of 2.09 Jy beam (the 1.4 mm data were 
re-imaged with the same beam-size as the 3.4 mm data), we get 
a spectral index of ~3.4. However, since our detection is barely 
above the 3<x level, and mml may well harbor a so far unde- 
tected hypercompact Hn region that could contribute still sig- 



nificant flux at 3.4mm (e.g.. IBeuther et al.ll2007al) . we refrain 
from further analysis of that feature. 

3.1 .2. The dense gas observed in CH 3 CN(5a: - 4 K ) 

In contrast to the 3.4 mm continuum emission, the CH3CN 
emission distribution shows the typical double-peaked 
morphology known from the previous NH 3 observations 
(IBeuther et alJl2005lEo07bh . The two CH 3 CN peaks appear to 
be associated with the two stro ngest mm continuu m sources 
(mml and mm2) in the region dHunter et al. I l2006h . Figure [3] 
presents the full CH3CN(5a: - 4#) spectra extracted toward the 
two peak positions, and clearly all five ^-components up to 
K = 4 with EJk = 128 K are well detected. Table [3] lists the 
fitted line parameters of the spectra. 
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50 
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50 




9.196 10 9.197 10 9.198 10 
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Fig. 3. CH 3 CN(5k - A K ) spectra (K = 0...4) extracted toward 
the two CH3CN peak positions in NGC6334I shown in Figure 
Q] The dotted line in the top panel shows a model spectrum 
created with XCLASS at a temperature of 200 K. 



Since the early work by lLoren & Mundyl (1 1984 CH 3 CN 
has often been used as a thermometer to estimate rotation 
temperatures of the dense gas via Boltzmann plots assuming 
optically thin emission in Local Thermodynamic Equilibrium 
(LTE). We tried this approach here as well, however, it 
failed because CH3CN(5a: - 4^) is optically thick. Similarly, 
we tried to model the CH3CN spectra in LTE using the 
XCLASS superset to the CLA SS software develop ed by Peter 
Schilke (priv. comm., see also IComito et al.ll2005l) . This soft- 
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Fig. 1. 3.4 mm spectral line and continuum emission toward NGC6334I: The grey-scale shows the 3.4 mm continuum and the 
full contours present the integrated CH3CN(54 - 44) emission from -12 to 2km s _1 . The dashed and dotted contours show the 
blue- and red-shifted HCN(l-O) emission with integration ranges [-25,-19] and [5,19]kms _1 , respectively. In the left panel, the 
3.4mm and CH3CN emission is contoured from 3cr and continue in 3<x steps (lcr(3.4mm) ~ 23 mJy beam -1 and lo^CHiCN) ~ 
7 mJy beam 1 ). The red- and blue-shifted HCN contours start at 2<x and continue in l cr steps (lcr(HCN - red) ~ 19 mJy beam 1 
and lcr(HCN - blue) ~ 16 mJy beam 1 ). The mm sources from iHunter e t al. (2006) are marked by stars. For clarity, the right 
panel shows the same data but with less contours, and additional poitions from various othe r observation s are in cluded (iden- 
tified at the bottom-left) . The corresponding references are: SMA mm continuu m data from lHunter et al.1 d2006h. H?Q maser s 
from lForster & Caswell! d 19891). CH 3 OH cl ass 11 masers from lWalsh et alj dl998l) . OH masers from lBrooks & Whiteoakl d200ll) . 
MIR sources from De Buize r et alj (2002 )). The NH3(6,6)/(8,6)/(l 1,9) are not shown but are spatially associated with mm2 
( Beuther et al.ll2007b : Walsh et al.ll2007 ). The synthesized beam is shown at the bottom right of each panel. 



ware package uses the line catalogs from JPL and CDMS 
dPovnter & Pickett! 1 9851 iMiiller et alj200lh . However, this ap- 
proach failed as well. Fig. [3] shows a model spectrum at T — 
200 K overlaid on the CH3CN(5a: - 4k) toward the main mm 
peak mml. While the K — 2, 3 components still fit relatively 
well, the model does neither reproduce the K - 0,1 nor the 
K = 4 component. While this is partly again an opacity ef- 
fect, it also shows that LTE is not appropriate for a source 
like NGC6334I. As outlined below, different /f-levels exhibit 
different line-widths and hence do not trace the same gas 
components. Temperature gradients within the sources are im- 
printed in the spectra further complicating single-temperature 
fits. More sophisticated modeling of the CH3CN emission is 
warranted. Although a few radiative tran sfer calculations of 



CH3CN data exist (e.g., lOlmi et al]l!993l) . the sparsely avail 



able collisional transition rates make such calculations a diffi- 
cult task (the only partly available data are a small compilation 
by Pei & Zengl 1995blah . Furthermore, in such hot and dense 
regions radiative excitation starts to matter which is hard to ac- 
count for in any modeling approach. 

While the Full Width Half Maximum (FWHM) line-widths 
Av for the two lowest energy lines of the K — 0, 1 compo- 
nents are less certain because of the line-blending between both 



components, one identifies a trend of increasing line-width, Av, 
with increasing K quantum number for K > 2 (Table [3). Figure 
|U shows the Av of CH 3 CN(5a: - A K ) for K > 2 toward the two 
CH3CN peaks in NGC6334I & I(N) plotted versus the level 
energy above ground, E u /k. For all four positions the trend of 
increasing Av versus E u /k is discernable. This trend indicates 
more internal motions from the warmer gas components (traced 
by the higher E u /k lines) which perhaps originate from the in- 
ner warm regions close to the protostars. Similar to that, the 
bottom-panel of Figure [5] shows the 2nd moment map of the 
CH3CN(54 - 44) line, i.e., its line-width distribution. Again we 
see the line-width increase toward the center close to the main 
mm continuum peaks. Different processes may cause such line 
broadening, e.g., accretion disk rotation, infall or outflow mo- 
tions. In particular, accretion disks are interesting candidates to 
explain such observational features. 

To investigate rotation from a potentially embedded mas- 
sive accretion disk, the top-panel of Figure [5] presents the 1st 
moment map of the CH3CN(54 - 44) line, i.e., its peak velocity 
distribution. Although this structure is only barely resolved by 
the synthesized beam of 2.5" x 1.8" we tentatively identify a ve- 
locity gradient across mml with an approximate position angle 
(PA) of 1 13 + 23 degrees from north. The same velocity struc- 
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Fig. 2. 3.4 mm spectral line and continuum emission toward NGC6334I(N) The four panels show the continuum (top-left), the 
integrated CH 3 CN(5 K ~^k) (K = 0, 1, top-right), and the red- and blue-shifted HCN( 1-0) emission (bottom-left and bottom-right 
respectively). The velocity regimes are given in the figure. The 3.4mm and CH3CN emission is contoured from 2<x and continue 
in 2<r steps (lcr(3.4mm) ~ 2.6 mJy beam -1 and lcr(CH3CN) ~ 4.9 mJy beam 1 ). The red- and blue-shifted HCN wing contours 
start at 2cr and continue in l<x steps (lcr(HCN - red) ~ 19mJybeam _1 and l<x(HCN - blue) ~ 15 mJybearrT 1 ). Negative 
features caused by insufficient uv-coverage are shown in dashed contours with t he same levels as th e emission. Markers of 
various other observations are identified in the top-right panel (SM A mm data from Hunter et al. (2006), the sources are labeled 
in the top-left panel, CH3OH class n maser from l Walsh et alj (1998), cm emission from lCarral et al.l(l2002h ). The two lines outline 



the axes of the two molecular outflows identified bv lMegeafh & Tieftrunki d 19991) Hunter et al. (in prep.). The synthesized beams 
are shown at the bottom right of each panel. 



ture is descernable in the lower K — 3,2 CH3CN transitions 
(we refrained from analyzing the K = 0, 1 lines because of 
their line-blending). To investigate this potential velocity gra- 
dient in more detail, we fitted the peak positions of each inde- 
pendent spectral channel (Fig. [6]) which should allow us to in- 
crease the resolving power to approximately 0.5 HPBW/(S/N), 
where HPBW equals the synthe sized beam and S/N the signal- 
to-noise ratio (Rei d et alJI 1988H . Similar to the moment map, 
the case is not clear-cut, but nevertheless, the data are indica- 



tive of a velocity gradient with a PA of ~ 134^ degrees 
from north. In comparison to these position angles, the PA of 
the previousl y identified molecular outflo w is ~ 46 degrees 
from north (Bachill er & Cernicharo 1990), which is approxi- 
mately perpendicular to that found from our measurements, es- 
pecially those from the highest-spatial-resolution position fit- 
ting (Fig.©. 

Based on these findings, we went back to the previous 
NH 3 (1,1) to (6,6) observations dBeuther et all |2005[ |2007b) 
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Table 3. Fitted CH 3 CN(5a; - 4 K ) line parameters 
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[km/s] 
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CH 3 CN(5 - 4 ) fl 
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CH 3 CN(5 3 - 4 3 ) 
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CH 3 CN(5 4 - 4 4 ) 
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NGC6334I, mm2 


CH 3 CN(5 - 4 )« 


25.6 


5.1±0.3 


CH 3 CN(5i -4i) n 
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" The fits to the K = 0, 1 lines are less accurate because of the strong 
line blending of both components. 
b Full Width Half Maximum (FWHM) 
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Fig. 4. Line-widths Av(CH 3 CN) from the various K- 
components plotted against the upper energy level of each line. 
We omitted the K = 0, 1 lines because of the line blending. The 
different symbols correspond to 4 different peak positions in 
NGC6334I and I(N) as labeled in the plot. The corresponding 
errors are given in Table [3] 

Adopting the proposed disk scenario, we can estimate 
the approximate rotationally supported binding mass M 10t as- 
suming equilibrium between the centrifugal and gravitational 
forces at the outer radius of the disk. Then we get 



searching for similar signatures in these data. While the main 
hyperfine lines as well as the satellite lines of NH3(1,1) are 
dominated by the large-sca le velocity gra dient over the two 



main cores (e.g., Fig. 13 in Beu ther et al.l 120051) . the satellite 



lines of the (J,K) lines with J,K>3 are overlapping and difficult 
to image. However, the satellite hyperfine lines of the NH 3 (2,2) 
transition exhibit, in addition to the velocity gradient over the 
two cores, a second velocity gradient across mml, again ap- 
proximately perpendicular to the large-scale outflow (Fig. |7}. 
The confirmation of this velocity gradient, first identified in the 
CH 3 CN(54 - 44) line, now in the lower opacity satellite hyper- 
fine NH 3 (2,2) line supports its general credibility. 

Since we do not resolve well the substructure of that veloc- 
ity gradient to better investigate the kinematics (e.g., is there 
any Keplerian motion present?) the data do not allow the claim 
that a massive accretion disk was detected. For example, an 
unresolved double-source with different lin e-of-sight velocitie s 



could produce a similar signature (e.g.. iBrogan et al.l 120071) . 
Nevertheless, these observations are suggestive of a rotating 
structure perpendicular to the molecular outflow within a pro- 
jected diameter of ~ 0.33" derived from Fig. [6] At the given 
distance of 1 .7 kpc, this corresponds to a radius of the rotating 
structure of -280 AU. This scale fits well to the sizes of the 
accretion disks simulated recently via 3-dimensional radiative - 
transfer hydrodynamic calculations bv lKnimholz et al. (2007). 



5v 2 r 
~G~ 



M rot [M ] = 1.13 10~ 3 Xc5v 2 [km/s] xr[AU] 



(1) 

(2) 



r is the disk radius, and 6v the Half Width Zero Intensity 
(HWZI) of the spectral line, approximately 5.1 km s _1 (half the 
velocity range shown in Fig. |6j. Equations Q] & [2] have to be 
divided by sin 2 (/) where i is the unknown inclination angle be- 
tween the disk plane and the plane of the sky (i = 90° for an 
edge-on system). With the given values we can estimate M rot to 
~ 8/(sin 2 (/)) M . This is of the sa me order as the mas s derived 
from the mm continuum emission (Hun ter et all2006b which is 
assumed to stem largely from the disk/envelope system. 

How do these masses correspond to the mass and lumi- 
nosity of the central embedded source? While the bolometric 
luminosity of NGC6334I is estimated to be < 2.6 x 10 5 L Q 
dSandel]||2o6ob . approximately 0.32 x 10 s L are attributed to 
the UC Hn region (based on the Lyman continuum flux pre- 
sented in IdePree et alJl995l) . One should keep in mind that this 
value may be a lower limit si nce dust could ab sorb a significant 
fraction of the uv-photons ( Kurtz et al. 19941) . The remaining 
< 2.3 x 10 5 L has to be due to the various sources associated 
with the hot molecular core. Since the associated mid-infrared 
source to the west (Fig. [Hi is of relatively low luminosity (only 
67 Ln. iDe Buizer et al.l 120021) . it is likely that most of the lu- 



minosity stems from the two main continuum and spectral line 
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Fig. 5. Moment maps of CH 3 CN(5 4 - 4 4 ) toward NGC6334I. 
The top panel shows the 1st moment (peak velocities) and the 
bottom panel the 2nd moment (line widths). The stars mark 
the positions of th e two main mm continuum sources from 
Hunte r et al. (2006). In addition to the velocity gradient from 
between mml and mm2, the 1st moment map exhibits a 2nd 
velocity gradient around mml in north-west south-east direc- 
tion (approximately perpendicular to the large-scale outflow). 



peaks. In the extreme case, splitting the luminosity simply by 
two, it still implies that about < 1.15 x 10 5 L o emanate from 
the mml region. Assuming a ZAMS star, this corresponds to 
an embedded star of <3OM . In the above adopted accretion 
disk scenario, this would imply an inclination angle i between 
the disk plane and the plane of the sky of approximately 30 
degrees. However, with the large uncertainties in the rotational 
mass and the central object mass estimate, such an inclination 
angle estimate should not be taken at face value, but only gives 
a rough idea that the system is neither edge- nor face-on. The 
comparable mass d erived from the mm continuum emission 
(Hun ter et al.l l2006) indicates that a significant fraction of the 
total system mass stems from the proposed accretion disk and 
envelope. This implies that the rotating structure is unlikely in 
Keplerian motion but that it may be a self-gravitating structure 
and potential site of ongoing sub-fragmentation (see also com- 
parable analytic calcula tion s and hydro-simulation s by, e.g., 
iKratter & Matznerll2006l and lKrumholz et al.1 d2007l) ). 

However, on a cautionary note, one has to keep in mind 
that the data are not conclusive as to whether we really see disk 
signatures or whether the observed gradient may be caused by 
other motions, e.g., an unresolved double-source. 
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Fig. 6. Positional offsets around mml of the different CH3CN 
(54-44) velocity channels derived via Gaussian fits of the peak 
emission in each separate channel. The error-bars are the nom- 
inal lcr errors from the fits, and the numbers label the central 
velocity for each position. 



3.1 .3. The molecular outflow observed in HCN(1-0) 

The blue- and red-shifted HCN(l-O) emission in Figure 
Q] shows high-velocity emission associated with the 
molecular outflow in north-east south-wester n direction 
dBachiller & Cernicharo|[l990l iLeurini et alj|2006l) . The PA of 
the emission is not exactly the 45 degrees derived previously 
from the single-dish CO observations but it is closer to 65 
degrees. However, such a discrepancy is not necessarily a sur- 
prise if one considers the missing flux problem we encounter 
in the HCN observations. Figure [8] (bottom panel) presents the 
HCN spectrum extracted toward the mml peak position, and 
while we see well the blue- and red-shifted emission, nearly 
all the flux around the Vi sr of ~ -7.6 km s _I is filtered out. This 
is also the reason why we do not see HCN emission from the 
hot core itself, which is prominent in CH3CN. Therefore, we 
just see some selectively chosen part of the outflow in HCN 
that could for example be associated with the limb-brightened 
cavity walls of the outflow which would explain the different 
apparent PA of the image. Figure [9] shows a position-velocity 
diagram centered at the main mm emission and CH3CN peak 
mml along the apparent axis of the HCN emission. Again we 
find no emission around the systemic vi sr but going to higher 
velocities, the blue- and red-shifted gas exhibits increasing 
velocity with increasing distance from the center, closely 
resem bling the typic al Hubble-law of molecular outflows 
(e.g., Lee et al. 200lb . Such Hubble-law like behavior can 
be explained on the smallest jet-scales close to the protostar 
by the decreasing gravitational potential of the central star, 
however, on the larger scales we observe here this effect 
gets negligible and the Hubble-law of molecular outflows 
is explained by a density gradient decreasing with distance 
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Fig. 7. First moment map (top) and position velocity diagram 
(bottom) of the the most blue-shifted satellite hyperfine line of 
NH?(2,2) toward N GC6334I. The data are re-examined from 
iBeuther et al.1 d2005l) . The position velocity digram is centered 
on mml with a position angle of 136 degrees from north. In 
addition to the velocity gradients between mml and mm2, the 
data exhibit a 2nd velocity gradient around mml in north-west 
south-east direction (approximately perpendicular to the large- 
scale outflow). The synthesized beam is 2.8" x 2.2". 



combined with the continuous (or episodic) driving of the jet 



that constantly induces energy in the outflow (e.g., IShu et al 
199lHSmith etai]ll997l:lDownes & Ravlll999h . 



Figure [10] (top panels) presents the 1st and 2nd moment 
map (intensity-weighted peak velocity and line-width distribu- 
tions) of the HCN(l-O) emission. It is interesting to note that 
the line-width distribution clearly peaks toward the strongest 
mm continuum and molecular line source mml and that the 
large-scale HCN velocity gradient is approximately centered 
toward that position as well. This is indicative of a scenario 
which puts the driving source of the molecular outflow ob- 
served in HCN at this position. However, previous NH?(6,6) 
maser observations showed the maser peak position being asso- 
ciated with mm2 with additional features along the axis of the 




NGC6334I 
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Av = 15.4 km s~' 
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main mm and CH 3 CN peak 
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Fig. 8. HCN(l-O) spectra in NGC6334I. The top spectrum 
shows the strong absorption toward the 3.4mm continuum 
peak that coincides with the UCHn region previously observed 
at cm wavelengths. The bottom spectrum is extracted toward 
the main CH3CN peak that coincides with the main mm con- 
tinuum peak SMA1 bv lHunter et al. I d2006h . 



CO molecular outflow (IBeuther et al. 2007b). These features 
were interpreted as indicative of this outflow being driven by 
a source associated with mm2. While the case is not clear-cut 
with two independent and different outflow driver indications, 
it is also possible that we are witnessing two different outflows 
observed in HCN and CO that may emanate from mml and 
mm2, respectively. In this scenario, mml could be the driver 
of a potentially denser and younger outflow that is better de- 
tected in HCN, whereas mm2 may be the driver of the larger, 
possibly older outflow observed in CO. Future observations at 
higher angular resolution and/or with different outflow tracers 
are required to assess the validity of this scenario. 

3.1.4. Absorption toward the UCHn region 

Another interesting feature of the HCN data is that we see 
strong absorption in the direction of the UCHn region, simi- 
lar to the previously observed absorption in CH3OH and NH3 
IBeuther et alJl2005l) . Figure [8] (top panel) shows the HCN(1- 
0) spectrum extracted toward the 3.4mm continuum peak po- 
sition. Fitting the HCN(l-O) line-width, we take into account 
its hyperfine structure consisting of three lines (F = — 1, 2 - 
1,1 - 1) with relative intensities of 1:5:3 in the optically thin 
limit and velocity shifts o f -7.1, and 4.9 km s -1 , respectively 
dPoynter & Pickett! fl985l) . The line- width Av of that absorp- 
tion feature is 15.4kms _1 , much broader than the previously 
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Fig. 10. Moment maps of HCN(l-O) toward NGC6334I. The left and right columns present the first and second moments (peak 
velocities and line widths), respectively. The top-row shows the moments for the emission of the outflow and ambient gas, whereas 
the bottom-row presents the corresponding parameters for the absorption features toward the UCHu region. The contours present 
the 3.4 mm continuum emission in 3cr contour levels (3cr ~ 69mJy beam -1 ). The markers are the same as in Fig.Q] 



observed line-widths between 1.7 and 2.1kms for CH3OH 
and between 1.3 and 1.9 kms" 1 for NH 3 (1,1) and NH 3 (2,2) 



Table 4. Line parameters for all absorption lines 



(Beu ther et al.ll2005l) . Table [4] lists the observed line-widths as 
well as the upper level excitation temperatures E u /k and the 
critical densities n CT i t for all observed absorption lines. The 
other extreme of the line-widths distribution is the line-width 
of the ionized gas o bserved in the H16a line Av = 32.0km s~' 
dde Pree et alJI 1995b . 



Inspecting the values in Table|4]one can discern two trends: 
the first is the distinction between the broad line- width and high 
critical density line HCN versus the small line-width and low 
critical density lines from NH3 and CH3OH. In addition, within 
the low critical density molecules, one can tentatively identify a 
correlation between increasing E lt /k and increasing line-width. 
Although judging from the errors this correlation is less clear, 
it remains suggestive. 



Line 




EJk 




Av 






[K] 


[cm- 3 ] 


[km/s] 


NH 3 (1,1) 




24 


2.1e3 


1.4 + 0.4 


NH 3 (2,2) 




65 


2.2e3 


2.0 + 0.3 


CH 3 OH(3 2 ,i 


-3,, 2 ) 


36 


1.0e4 


1.8 ± 0.1 


CH 3 OH(4 2 , 2 


-4,, 3 ) 


45 


1.2e4 


2.1 +0.1 


CH 3 OH(2 2 , 


-2 U ) 


29 


0.8e4 


1.7 ±0.1 


HCN(l-O) 




4 


2.6e6 


15.4 ±0.3 



The NH 3 and CH 3 OH data are from Beut her et alj d2005l) . The listed 
parameters are the upper level excitation energies E u /k, the critical 
densities n cl jt calculated at 60 K and the observed line-widths of the 
absorption lines toward the UCHu region peak position. 



There are different possibilities to e xplain such trends: T he 
picture of an expanding UCHu regiorQ ( Beuther et al. 2005) in 



1 In contrast to the previous CH 3 OH and NH 3 absorption lines that 
showed additional blue-shifted emission indicative of expanding gas 
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Fig. 9. Position- velocity diagram of HCN(l-O) in NGC6334I 
along the main outflow axis. The cut is centered on the main 
CH 3 CN peak at positional offset -0.3"/4.5" with a position 
angle of 65° from north-to-east. Emission at the velocity of rest 
around -7.6kms _I is largely filtered out. 



its surrounding envelope implies that closer to the UCHn region 
surface the molecular densities and temperatures are higher 
than further outside. Therefore, spectral lines tracing higher 
densities around an expanding UCHn region are expected to 
exhibit broader line-widths because their emitting gas is more 
directly impacted than the lower-density medium further out. A 
similar explanation could hold for the different excitation tem- 
perature regimes as well. A point of caution is that the opacity 
of HCN is probably one to two orders of magnitude larger than 
that of, e.g., NH3. While we do not detect any NH3 satellite hy- 
perfine structure lines in absorption implying low NH3 opaci- 
ties, we cannot infer that exactly for HCN. If the HCN(l-O) op- 
tical depth were that high that it could not trace the dense gas 
regions close to the expanding UCHn region, then the above 
picture could hardly hold. Another explanation is based on 
the enhancement of HCN in the molecular outflow. Although 
the outflow does not emanate from the U CHn region but from 
the neighboring mm continuum sourcefs Ude Pree et al.l (119951) 
found a velocity gradient in the ionized gas similar to that of 
the molecular outflow. They suggest that the molecular out- 
flow^) in the region may well disturb the velocity field of the 
UCHn region and produce the velocity gradient this way. Since 
HCN is known to b e abundant in mol ecular outflows (see also 
IRAS 18566+0408. IZhang et al.ll2007[ and IRAS 20126+4104, 
Liu et al. in prep.) a similar effect could explain the broad 
HCN(l-O) line-widths toward the UCHn region. However, the 
1st and 2nd moment maps (peak velocities and line- widths) of 
the HCN(l-O) absorption presented in Fig. [TOl (bottom panels) 
show no velocity gradient but a peak of the line-width distri- 
bution close to the center of the UCHn region. This is counter- 
intuitive for the outflow picture and suggests that the line-width 



around the UCHn region (Beu ther et alj|2005l) . the HCN data do not 
exhibit any clear expansion or infall signature. 



differences are caused by the expanding UCHn region, possibly 
in a fashion comparable to that proposed above. Nevertheless, 
without knowing the optical depth of the HCN(l-O) line we 
cannot distinguish between the two scenarios. 

3.2. NGC6334I(N) 

3.2.1. Millimeter continuum emission 

Figure [2] (top-left panel) presents the 3.4mm continuum emis- 
sion toward NGC6334I(N). In contrast to NGC6334I, where 
we only see the UCHn region in the 3.4 mm continuum, to- 
ward NGC6334I(N) we detect four out of seven previously 
identified protostel lar 1.4 mm dust continuum condensations 
( Hunter et al. 20061) above a 5cr level of 13 mJybeairT 1 . The 
peak fluxes 5 pea k(3.4mm) of the four sources, mml to mm3 
and mm6, are listed in Table [5] Emission features below 5<x 
are not considered further. The 3.4mm peak associated with 
mml is the strongest and shows an additional extension to- 
ward the north-western 1 .4 mm source mm5. In contrast to that, 
the 1.4 mm peak mm4, which is associated with cm continuum 
and CH3OH class 11 maser emission, is not detected in our data 
above the 5<x level. 

For a better comparison of the 3.4 mm fluxes with the previ- 
ously observed 1 .4 mm observations, we re-imaged the 1 .4 mm 
data with exactly the same synthesized beam of 2.6" x 1.8" (po- 
sition angle of 84 degrees from north). Figure [TT] presents an 
overlay of the SMA 1.4mm data with the ATCA 3.4mm data 
at this same spatial resolution. While mml and mm6 are clearly 
separated in both wavelength bands, it is interesting to note that 
mm2 and mm3 merge in the 1 .4 mm image at the reduced lower 
spatial resolution. The corresponding 1 .4 mm peak fluxes are 
listed in Table[5] For the well separated sources mml and mm6 
we can estimate the spectral indices a between both bands now, 
the derived v alues are 3.7 and 3.1, re spectively (Table[5]i. 

Recently Rodriguez et al. ( 20071) showed that the cm emis- 
sion from NGC6334I(N) is caused by free-free emission, 
whereas shortward of 7 mm wavelength the spectral energy 
distribution is dominated by dust continuum emission. In the 
Rayleigh- Jeans limit, the flux S scales with S oc y 2+ ^ where ft 
is the dust opacity index. With the measured spectral index a, 
we have dust opacity indices of 1 .7 and 1 . 1 for mml and mm6, 
respectively. While the (3 value of mml is consistent with of- 
ten observed values between 1.5 and 2, it is lower for mm6. 
Although the synthesized beams of both maps are the same, 
the uv-coverage was not during the observations. Furthermore, 
the continuum maps show a more peaked morphology toward 
mml than toward mm6. Thus, it is possible that mm6 has 
more extended structure and that this may be sampled better by 
the ATCA observations, possibly accounting for the observed 
lower values of a and /3 toward mm6. However, it is also feasi- 
ble that physical reasons are responsible for this effect because 
a decreasing j3 can be attributed to grain growth in circumstellar 
disks as well (e.g. jBeckwith et al.ll 199QD . Nevertheless, in this 
scenario, it appears surprising that mml, which likely drives 
an outflow and hence probably contains an accretion disk, has 
a value of j3 close to the interstellar medium values, whereas 
the source mm6, that does not exhibit clear signs of ongo- 
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ing star formation, shows a lower value. Low a values could 
in principle also be caused by high optical depth or an early 
high-frequency turnover of the Planck-function caused by the 
low temperatures. However, the Planck-turnover of even a cold 
core starts changing the slope of the Rayleigh-Jeans part of the 
spectral energy dis tribution usually only above 350 GHz (e.g., 
iBeuther et al.ll2007ah . and high optical depth compared to mml 
appears unlikely as well. Therefore, spatial filtering as well as 
physical effects could produce the low a toward mm6, but we 
cannot set better constraints here. 

It should be noted that the spectral indices we find be- 
tw een 3.4 and 1.3 mm are larger than those derived previously 
bv iRodrfguez et al. (2007) ba sed on VLA 7 mm data and the 
1 .3 mm data by iHunter et al.1 d2006l) . They find a of -2.4 for 
both sources mml and mm6. A potential explanation for this 
discrepancy can arise from different dust components traced 
by the various arrays. While the VLA 7 mm observations could 
be dominated by a compact dusty disk-component that remains 
optically thick at shorter wavelengths, the 3.4 and 1.4mm data 
are likely dominated by the larger-scale optically thin enve- 
lope. This envelope with spectral indices between 3.1 and 3.7 
(for mml and mm6) is weak and below the detection limit at 
7 mm, which can cause the non-detection of any large-scale 
emission in the VLA data. A general word of caution should 
be added that comparing fluxes from different interferometers 
may be unreliable because spatial filtering usually affects the 
fluxes of most measurements. Reducing datasets with similar 
uv-coverage and the same synthesized beam, as we have done 
here for the 3.4 and 1.4mm data, helps to minimize this prob- 
lem. 

Table 5. Millimeter continuum parameters for NGC6334I(N) 



Source 


5 pea k(3.4mm) 


Speak(1.4mm) 


a 




[mjy/beam] 


[mjy/beam] 




mml 


36.6 


1077 


3.7 


mm2 


14.2 


524° 




mm3 


14.2 


524" 




mm6 


17.1 


285 


3.1 



' The sources mm2 and mm3 are merged in the 1 .4 mm dataset. 



3.2.2. Dense gas traced by CH 3 CN 

The dense gas traced via the CH3CN(5# - 4^) /f-ladder ex- 
hibits two prominent peaks associated with the two main mm 
continuum sources mml and mm2 (Fig. [2J. There is an ad- 
ditional tentative 3<x CH3CN peak toward the mm peak mm4 
associated with the cm and CH3OH maser position, however, 
we refrain from further interpretation of this because the map 
shows negative features due to the incomplete uv-sampling and 
hence poor deconvolution on a comparable level. The other 1 .4 
and 3.4 mm continuum sources are not detected in the CH3CN 
emission. This indicates that the average densities and/or abun- 
dances at those other mm positions are likely lower than the 
critical densities of the CH3CN(5# - 4k) lines of a few times 
10 s cirT 3 (TableQ}. Based on this, one can speculate that except 
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Fig. 11. Overlay of the SMA 1.4 mm continuum map (grey- 
scale with dashed contours) with the ATCA 3.4mm map (solid 
contours) at the same angular resolution of 2.6" x 1.8". The 
1 .4 mm data are contoured from 10 to 90% of the peak emission 
(1077 mJy beam -1 ). The 3.4mm emission is contoured in 2cr 
steps (2cr(3.4mm) ~ 5.2 mJy beam -1 . The synthesized beam of 
both images is the same with 2.6" x 1.8" (position angle of 84 
degrees from north). 



for the sources mml, mm2 and mm4, all other mm continuum 
peaks may still be in a pre-stellar phase prior to active star for- 
mation. 

Figure Q~2] presents the whole CH3CN /f-ladder spectrum 
extracted from the data-cube toward the two peak positions, 
and we again detect all ^-components in both sources. The 
spectrum toward mml suffers from very high opacity prob- 
lems as already outlined for NGC6334I ( CTQl ), prohibit- 
ing any temperature estimate. However, for mm2 the situa- 
tion is slightly different and temperature estimates become 
feasible. Again using the XCLASS software, we produced 
model spectra, and the CH3CN(5a: - spectrum can be fit- 
ted with temperatures in the regime of 170 + 50 K. Such high 
temperatures additionally confirm the hot core-like nature of 
these sub-sources w ithin the over all still relatively cold region 
NGC6334I(N) (e.g.. lGezari|[l982l) . 

Since previous NH;(6,6) obs ervations showed a double- 
horn spectral profile toward mml (IBeuther et al . 2007b), which 
could potentially be produced by an underlying accretion disk, 
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Fig. 12. CH 3 CN(5 / r : - A K ) spectra (K = 0...4) extracted to- 
ward the two CH 3 CN peak positions in NGC6334I(N) shown 
in Figure [2] The dotted line in the lower panel shows a model 
spectrum created with XCLASS at a temperature of 170 K. 



we searched the high-Zf CH3CN lines for rotation signatures as 
well. Figure [13] presents the 1st and 2nd moment maps of the 
CH3CN(54 - 44) line toward mml. While we see a line-width 
increase toward the center, which is indicative of more activity 
in the central region (e.g., rotation, outflow, infall), we cannot 
identify any obvious velocity gradient. Therefore, the data do 
not allow a more detailed analysis of the previously proposed 
underlying disk. 

3.2.3. A processing outflow? 

The HCN(l-O) emission suffers strongly from the missing 
short spacings as indicated by all the negative features seen 
in Figure [2] The morphology of the remaining HCN emis- 
sion structures nevertheless indicates that most of it is asso- 
ciated with the north-east south-west molecular outflow than 
with the perpendicular one. Figure fT4l shows the 1st and 2nd 
moment maps of the HCN emission, and the peak velocity 
and line-width distributions are a bit peculiar. The most red- 
and blue-shifted emission is centered around mml indicating 
that this is likely the driver of the north-east/south-west out- 
flow. However, going further to the north-east along the ma- 
jor outflow axis we find an inversion of the outflow velocities 
from red-shifted close to the protostar to blue-shifted at offsets 
around (12"/10") from the phase center and then again red- 



o 
o 
Q 

< 



5 - 



4 - 



3 - 



CH 3 CN(5 4 -4 4 ) 



|[km/ s] 




loment 1 : Peak velocities 



I I | I I I I | I I I I | I I I 
CH 3 CN(5 4 -4 4 ) 




■ 



[km/ s] 

6 



loment 2: Line widths 
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Fig. 13. Moment maps of CH 3 CN(5# - A K ) toward 
NGC6334I(N) mml. The top panel shows the 1st mo- 
ment (peak velocities) and the bottom panel the 2nd moment 
(line widths). The star marks the positions of the main mm 
continuum sources from Hunt er et al.1 (|2006). 



shifted emission further outward at offsets around (17"/ 13"). 
Regarding the line-width distribution, the 2nd moment map in 
Figure [14] shows that the red-shifted parts of the north-eastern 
outflow wing have systematically less broad lines than the blue- 
shifted part. This whole behavior is less pronounced in the 
south-western outflow direction. 

While the line-width differences are hard to explain, the 
change between blue-and red-shifted outflow emission is con- 
sistent with a picture in which a mean molecular outflow axis is 
close to the plane of the sky, and where the outflow precesses 
around that axis, producing in some region the receding red- 
shifted and in other regions the approaching blue-shifted emis- 
sion. This precession picture is further supported by the spatial 
distribution of the HCN emission which resembles a reverse 
S-shaped morphology in the 1st moment map (Fig.[T4li. 

The 2nd south-east north-west outfl ow is so far only ob- 
served at lower spatial resolution by iMegeath & Tieftrunk 
who did not identify a driving source for it. While 
it could emanate from any other source, e.g. mm2, it is also 
possible that it originates from the vicinity of mml. In the 
framework of precessing jets, a likely reason to cause the pre- 
cession is the existence of multiple embedded objects (e.g., 
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Fendt & Zinnecker 1998) which hence could also produce the 



quadrupolar outflow morphology (e.g. JGueth et al 



2001). 
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Fig. 14. Moment maps of HCN(l-O) toward NGC6334I(N). 
The top panel shows the 1st moment (peak velocities) and the 
bottom panel the 2nd moment (line widths). The symbols cor- 
respond to the same phenomena as in Fig. [2] 



4. Conclusions and Summary 

Our 3.4mm continuum and CH3CN/HCN spectral line study 
of the massive twin cores NGC6334I and I(N) reveals many 
new insights into that intriguing pair of massive star-forming 
regions. Both sets of spectral lines as well as the continuum 
emission are clearly detected toward both targets. While the 
continuum emission in NGC6334I mainly follows the UCH11 
regions and the strongest protostellar 1.4 mm peak is detected 
only at a ~ 3<x level, in NGC6334I(N), the 3.4mm continuum 
emission traces four of the previously identified protostellar or 
pre-stellar condensations. 

In both regions, the whole CH3CN(5a: — 4k) A'-ladder from 
K = to 4 is detected toward the strongest protostellar conden- 
sations. While the emission is in most cases so optically thick 
that temperature estimates are prohibited, toward the secondary 
mm peak in NGC6334I(N) we can estimate a temperature of 
170 + 50 K. Toward all four detected CH3CN emission sources, 



we find a correlation between increasing line-width and in- 
creasing excitation temperature of the K components. Since in- 
creasing excitation temperatures are expected closer to the pro- 
tostars, this implies more internal motions, e.g., outflow, infall 
or rotation, the closer one gets to the central protostar. Similar 
signatures are observed in the CH3CN 2nd moment maps. 

To investigate potential rotation, we produced 1st moment 
maps toward all CH3CN peaks, and fitted the channel peak po- 
sitions toward mml in NGC6334I to effectively increase the 
spatial resolution. We identify a velocity gradient toward mml 
in NGC6334I that is oriented approximately perpendicular to 
the known large-scale outflow. This may be interpreted as a 
signature of a rotating structure, maybe associated with a mas- 
sive accretion disk. While early rotation-disk c laims for that 



regio n were on scales of the molecular outflow dJackson et al 
Il988l) . we are now reaching spatial scales of the order a few 
hundred AU, much more reasonable for accreti on disks (e.g., 
lYorke & SonnhalterteOollKrumhoIz et al.l2007l) . However, we 
have to stress that the rotation signatures are not conclusive yet 
because, e.g., an unresolved double-source could produce simi- 
lar signatures. Further investigation at higher angular resolution 
are required to resolve this issue. While higher angular resolu- 
tion images with the ATCA at 3 mm wavelengths are difficult 
because of decreasing phase stability with increasing baseline 
length, one may tackle that problem in some highly excited 
NH3 lines in the 12 mm band. Furthermore, ALMA will allow 
the investigation of this source in much greater depth. Toward 
the previously suggested disk candidate in NGC6334I(N) we 
cannot identify a similar velocity gradient. 

In contrast to conventional wisdom that HCN traces the 
dense gas cores, we find it most prominently in the molecu- 
lar outflows of both massive star-forming regions. The velocity 
structure of the outflow in NGC6334I is relatively normal and 
follows the well-known Hubble-law for molecular outflows. In 
addition to that, we find the broadest HCN line-width toward 
the main mm continuum peak mml. In contrast to the previ- 
ously found elongation of NH3(6,6) maser emission indicating 
that mm2 could drive the molecular outflow, these data sug- 
gest that mml harbors the outflow driving source. However, 
it is also possible that there are two molecular outflows with 
different properties that are preferentially detected in different 
tracers (HCN in this work and CO for the previous larger-scale 
outflow detection at a slightly different position-angle). The ve- 
locity structure of the NGC6334I(N) outflow is more peculiar. 
There we find a change between blue- and red-shifted outflow 
emission on one side of the outflow. Taking into account the ad- 
ditionally bended morphology of that outflow, a possible expla- 
nation for this velocity structure is a precessing outflow close 
to the plane of the sky. 

Furthermore, HCN exhibits a broad absorption feature with 
a line-width of ~ 15.4km s 1 toward the UCH11 region in 
NGC6334I. Comparing the line-width of the previously ob- 
served absorption features of NH3 and CH3OH, which are of 
the order 2 km s _1 , with the HCN line-width as well as the line- 
width observed in the ionized gas of ~32kms _1 , two expla- 
nations are possible to explain the different line-widths. The 
velocity gradient identified in the ionized gas indicates that 
it may be influenced by the molecular outflow close by. If 
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the opacity of HCN were very large it traced only the outer 
gas layers around the UCHn region and could also be influ- 
enced by the outflow. On the other hand, NH3 is optically thin 
based on the absent absorption in the hyperfine satellite lines. If 
HCN were optically thin as well, another possible explanation 
would be based on the different critical densities of the var- 
ious molecules: In the picture of an expanding UCHn region, 
the densities close to the UCHn region surface should be higher 
than those further outside. Hence HCN may trace the gas closer 
to the expanding UCHn region and is then much stronger af- 
fected by the expansion process that the lower-density regions 
further out that are traced by NH3 and CH3OH. To differenti- 
ate between both models, observations of a rarer HCN isotopo- 
logues are required to derive its optical depth. 
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